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^«0 ■ Abstract. We discuss the role of turbulence in establishing the observed emission 

measures and ionization of the warm ionized medium. A Monte Carlo radiative transfer 
code applied to a snapshot of a simulation of a supernova-driven, multi-temperature, 
stratified ISM reproduces the essential observed features of the WIM. 
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In recent years, it has become clear that turbulence plays a central role in estab- 
lishing the observed conditions in the warm ionized medium (WIM) of the Galaxy. 
Although the WIM was first de tected four decades ago, the 1996 completion of the 

O I Wisconsin Ha Mapper (WHAM: lHaffner et al.ll2003h made possible large-scale studies 

of the physical conditions in the medium, particularly through the combination of Ha 

^ I and coUisionally excited me tal lines. These and o ther studies have established that the 

WIM is near ly fully ionized dReynolds et al.ri998h . primarily by the O star Lyman con- 
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tinuum flux (lReynoldslll990l) . However, the mechanism by which ionizing photons es- 



^ ■ cape the immediate vicinity of the O st ars near the Galactic plan e to reach the observed 

r^ ; 1 .0 - 1 .5 kpc scale height of the WIM ("Savage & Wakkei^' 20091') has been the subject of 

^ ■ numerous investigations (Miller & Cox 1993; Dove & Shulll ll994 . ICiardi et alT 



20021 : 



Wood & Mathisll2004l : iHaffner et al.ll2009l) . These studies have generally found that 



the propagation of ionizing photons to \z\ ~ 1 kpc in an ISM with the observed densi- 
ties requires the presence of low-density paths through the ISM, but none produce the 
Q ' required low-density paths in a physically motivated way. 

Meanwhile, the WIM has most often been modeled as consisting of either pho- 
toionized shells surrounding neutral clouds (McKee & Ostriker 1977) or as discret e 



clouds occupying some fraction of the total volume of the Galaxy (lReynoldslll99lh . 
^ I However, this picture is clearly incomplete. The WIM has a hig h Reynold s number and 

therefore is turbulent (Benjamin 1999 ). as seen observationally (lArmstrong et al. 1995) . 
The recent advent of 3D simulations of supernova-driven tu rbulence in a multiphase 
ISM \yhich cover dyn amically significant time and size scales dAvillez & Breitschwerdtl 



20041 : iJoung & Mac Low..20 06) has made possible detailed studies of the interplay be- 



tween the physical conditions revealed by optical emission lines and turbulence, helping 
to address the long-standing puzzle of the ionization of the medium. 



1. Distribution of emission measure 

First, we consider the global distribution of Ha emission measures observed in the 
WHAM Northern Sky Survey dHaffner et al.ll2003l ). We derived the emission measure, 
EM = j nl ds, from the Ha surface brightness assuming a temperature of 8000 K. We 
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Figure 1 . Grey lines: Histogram of EM sin \b\ for diffuse WIM sightlines with 
\b\ > 10° from the WHAM survey, duplicated in all three panels. Black lines: His- 
togram of simulated EM derived from Ms ~ 0.7 (left), Ms ~ 2 {center), and Ms ~ 7 
(right) externally driven, isothermal MH D simulations. D ashed lines: lognormal fits 
to the simulated EM distributions. From lHill et alj (I2008I) . 



then chose a subset of the WHAM survey which excludes sightlines towards large H II 
regions associated with identified, discrete ionizing stars. We also excluded sightlines 
with Galactic latitude \b\ < 10°, where extinction due to dust may affect the measured 
Hot surface brigh tness. T he distribution of EM sin \b\ is well approximated by a lognor- 
mal distribution (IHiU et al...2008„ and Figure [T])- 



2. Simulations 



A random series of compressions and rarefactions, such as those that would occur in a 
turbulent fluid, establish a lognormal distribution of density (Vazquez-Semadeni 1994). 
This l ed us to consider 3 D MHD simulations of artificially driven, isothermal turbu- 
lence dKowal et al.ll2007h . These simulations are generic and dimensionless. We ap- 
plied physica l scales appropri ate to the WIM and derived simulated emission measure 
distributions dHill et al.ll2008L and Figure[T]). In all cases, the simulated distributions are 
approximately lognormal. The width of the distribution is highly sensitive to the sonic 
Mach number (M.,) of the turbulence and is therefore the best parameter to compare to 
the observations. Subsonic models (M^ ~ 0.7) produce a narrower EM distribution than 
observed, while highly supersonic (Mg « 7) models produce a wider EM distribution 
than observed. However, the mildly supersonic (M^ « 2) simulations produce an EM 
distribution consistent with observations. 

Although the isothermal, externally driven simulations repr oduce the observed 
EM d istrib ution in the WIM , the driving mechanism is artificial. iJoung & Mac Lovvl 
(l2006l) and lJoung et all (l2009h developed 3D hydrodynamical simulations of a 0.5 kpcx 
0.5 kpcx 10 kpc region of the ISM with parameterized heating and cooling and a Galac- 
tic gravitational potential. Supernova explosions stir the medium, driving turbulence 
and structure formation. These simulations produce a multi-phase, stratified ISM which 
is well suited for a more complete study of the role of turbulence in the WIM. 

To investigate the WIM in the context of a supernova-driven turbulent ISM, we 
use a 3D Monte Carlo photoionization code (Wood & Loeb 2000) to model the pho- 
toionization of gas in a snapshot of the Joung & Mac Low supernova-driven models. 
We randomly place six O stars in the simulated, initially neutral ISM with a Gaussian 



Turbulence in the WIM 




Figure 2. Mean density of neutral (solid grey line) and ionized (solid black line) 
gas in the hydrody namical-photoionization simulations as a function of height. 
iDickey & LockmanI (dashed grey line) and exponential (dash ed black line) de nsity 
distributions for neutral and ionized gas are also shown. From lWood et al.l (l2010i) . 



distribution and a scale height of 63 pc. The surfa ce density of 24 O sta rs kpc^ is cho- 
sen to represent a sample of O stars near the sun (JGarmany et al.|[l982l) . We distribute 
the flux evenly between the O stars, investigating total ionizing luminosities Q over the 
range 10"^*^ s"^ < 2 < 10^" s~^ This span reflects a variation in efficiency of ionizing 
photons escaping from the immediate v icinity of the s tars; the total ionizing flux 
available from six O stars is » 10^° s~' (|Vacca et alJll996h . While neutral hydrogen is 
opaque to Lyman continuum radiation, the optical depth of nearly fully ionized hydro- 
gen is much lower. Therefore, we use an iterative process to establish a sel f-consistent 
distrib ution of neutral and ionized hydrogen. Our results are presented by IWood et al.l 
(1201 Ol) and summarized here. For ionizing luminosities Q >2x 10^^ s"\ ionizing pho- 
tons reach the top of the grid at |z| = 2 kpc. Figure |2] shows the horizontally-averaged 
density of the neutral and ionized gas as a function of height for the Q = 2 x 10^^ s"^ 
case. At heights |z| > 1 kpc, the gas is nearly fully ionized. 



3. Conclusions 



Supernova-driven turbulent simulations combined with photoionization modeling pro- 
duces the most se lf-consist modeling of t he warm ionized med ium to date. This finding 
is consistent with lCiardi et al.l(l2002h and lHaffner et al.l(l2009h . who found that analytic 
fractal density structures produce sufficient low-density pathways for the ionizing pho- 
tons to travel large distances. In combination with energy budget arguments (Reynolds! 
19901) . these results strongly suggest that turbulence allows photoionization from hot 
stars to be the primary ionizing mechanism in the WIM, whereas a traditional picture 
of discrete clouds in the ISM does not allow the necessary propagation of ionizing 
photons at realistic ISM densities. 

Inconsistencies between the combined hydrodynamical and photionization simu 



lations and the observations suggest avenues for future investigation. The Joun g et al. 



(12009) hydrodynamical simulation produ ces a more centrally peake d gas distribution 
than observed, with densities higher than a lDickey & Locknianl(ll990 ) distrib ution near 
the midplane and lower than Dickey-Lockman at |z| ~ 0.2 kpc (Fig. [2] and iWood et al. 
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2010'), possibly partly due to effects of resolution or box size dAvillez & Breitschwerdj 



2004). The inclusion of magnetic fields in our supernova-driven simulations (work in 
progress) may help to address this inconsisten cy by applying additional vertical pres- 
sure, although lAvillez & Breitschwerdtl (l2005b found that the magnetic field in their 
simulations is largely uncorrected with density. 

Also, although our simulations reproduce the existence of the WIM at large heights 
without resorting to an ad hoc interstellar density structure, the distribution of emission 
measures is considerably wider than in the observations o r the mildl y supersonic artifi- 



cially driven, isothermal simulations shown in Fig. [Tland lHill et al.l (2008,). Therefore, 
the hydrodynamic supernova-driven simulations produce more contrast between low 
and high density gas than obse rved in the Galaxy. Because magnetic fields reduce the 
compressibility of the medium (l Avillez & Breitschwerdt,2005.) . we expect the inclusion 
of magnetic fields in the simulations to improve this inconsistency as well. 
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